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Abstract 

There has recently been considerable interest in the possibility, 
both theoretical and practical, of invisibility (or "cloaking" ) from ob- 
servation by electromagnetic (EM) waves. Here, we prove invisibility 
with respect to solutions of the Helmholtz and Maxwell's equations, 
for several constructions of cloaking devices. The basic idea, as in the 
papers |GLU21 IGLU31 [Lei IPSSl] , is to use a singular transformation 
that pushes isotropic electromagnetic parameters forward into singu- 
lar, anisotropic ones. We define the notion of finite energy solutions 
of the Helmholtz and Maxwell's equations for such singular electro- 
magnetic parameters, and study the behavior of the solutions on the 
entire domain, including the cloaked region and its boundary. We 
show that, neglecting dispersion, the construction of |GLU3[ IPSSl] 
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cloaks passive objects, i.e., those without internal currents, at all fre- 
quencies k. Due to the singularity of the metric, one needs to work 
with weak solutions. Analyzing the behavior of such solutions inside 
the cloaked region, we show that, depending on the chosen construc- 
tion, there appear new "hidden" boundary conditions at the surface 
separating the cloaked and uncloaked regions. We also consider the 
effect on invisibility of active devices inside the cloaked region, inter- 
preted as collections of sources and sinks or internal currents. When 
these conditions are overdetermined, as happens for Maxwell's equa- 
tions, generic internal currents prevent the existence of finite energy 
solutions and invisibility is compromised. 

We give two basic constructions for cloaking a region D contained 
in a domain 17 C M",n > 3, from detection by measurements made 
at dn of Cauchy data of waves on J7. These constructions, the single 
and double coatings, correspond to surrounding either just the outer 
boundary dD^ of the cloaked region, or both dD~^ and dD~, with 
metamaterials whose EM material parameters (index of refraction or 
electric permittivity and magnetic permeability) are conformal to a 
singular Riemannian metric on fi. For the single coating construction, 
invisibility holds for the Helmholtz equation, but fails for Maxwell's 
equations with generic internal currents. However, invisibility can be 
restored by modifying the single coating construction, by either in- 
serting a physical surface at dD~ or using the double coating. When 
cloaking an infinite cylinder, invisibility results for Maxwell's equa- 
tions are valid if the coating material is lined on dD~ with a surface 
satisfying the soft and hard surface (SHS) boundary condition, but in 
general not without such a lining, even for passive objects. 



1 Introduction 

There has recently been considerable interest |AEl IMNt [Let IPSSll IMBW] 
in the possibility, both theoretical and practical, of a region or object be- 
ing shielded (or "cloaked") from detection via electromagnetic (EM) waves. 
|GLUll §4] established a non-tunnelling result for time-independent Schrodinger 
operators with highly singular potentials. This can be interpreted as cloaking, 
at any frequency, with respect to solutions of the Helmholtz equation using a 
layer of isotropic, negative index of refraction material. |GLU2l[GLU3j raised 



the possibility of passive objects being undetectable, in the context of electri- 
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cal impedance tomography (EIT). Motivated by consideration of certain de- 
generating families of Riemannian metrics, families of singular conductivities, 
i.e., not bounded below or above, were given and rigorous results obtained for 
the conductivity equation of electrostatics, i.e., for waves of frequency zero. 
A related example of a complete but noncompact two-dimensional Rieman- 
nian manifold with boundary having the same Dirichlet-to-Neumann map as 
a compact one was given in [LTUj . 

More recently, there has been exciting work based on the availability of meta- 
materials which allow fairly arbitrary behavior of EM material parameters. 
The constructions in |Lej are based on conformal mapping in two dimensions 
and are justified via change of variables on the exterior of the cloaked region. 
[PSSlj also proposes a cloaking construction for Maxwell's equations based on 
a singular transformation of the original space, again observing that, outside 
the cloaked region, the solutions of the homogeneous Maxwell equations in 
the original space become solutions of the transformed equations. The trans- 
formation used is the same as used previously in |GLU2l IGLU3] in the con- 
text of Calderon's inverse conductivity problem. The paper |PSS2] contains 
analysis of cloaking on the level of ray-tracing, while full wave numerical sim- 
ulations are discussed in [CPSSPj . Striking positive experimental evidence 
for cloaking from microwaves has recently been reported in |SMJCP^ . 

Since the metamaterials proposed to physically implement these construc- 
tions need to be fabricated with a given wavelength, or narrow range of 
wavelengths, in mind, it is natural to consider this problem in the frequency 
domain. (As in the earlier works, dispersion, i.e., dependance of the EM 
material parameters on k, which is certainly present for metamaterials, is 
neglected.) 

The question we wish to consider is then whether, at some (or all) frequencies 
k, these constructions allow cloaking with respect to solutions of the Helm- 
holtz equation or time-harmonic solutions of Maxwell's equations. We prove 
that this indeed is the case for the constructions of [GLUSt IPSSlj . as long 
as the object being cloaked is passive; in fact, for the Helmholtz equation, 
the object can be an active device in the sense of having sources and sinks. 
On the other hand, for Maxwell's equations with generic internal currents, 
invisibility in a physically realistic sense seems highly problematic. We give 
several ways of augmenting or modifying the original construction so as to 
obtain invisibility for all internal currents and at all frequencies. 
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Due to the degeneracy of the equations at the surface of the cloaked re- 
gion, is important to rigorously consider weak solutions to the Helmholtz 
and Maxwell's equations on all of the domain, not just the exterior of the 
cloaked region. We analyze various constructions for cloaking from obser- 
vation on the level of physically meaningful EM waves, i.e., finite energy 
distributional solutions of the equations, showing that careful formulation of 
the problem is necessary both mathematically and for correct understanding 
of the physical phenomena. It turns out that the cloaking structure imposes 
hidden boundary conditions on such waves at the surface of the cloak. When 
these conditions are overdetermined, finite energy solutions typically do not 
exist. The time-domain physical interpretation of this is not entirely clear, 
but it seems to indicate an accumulation of energy or blow-up of the fields 
which would compromise the desired cloaking effect. 

As mentioned earlier, the example in |PSSlj turns out to be a special case 
of one of the constructions from |GLU2t IGLU3 J. which gave, in dimensions 
n > 3, counterexamples to uniqueness for Calderon's inverse problem [Cj for 
the conductivity equation. (Such counterexamples have now also been given 
for n = 2 [Vl IKSVW] .) Thus, since the equations of electromagnetism (EM) 
reduce at frequency = to the conductivity equation with conductivity 
parameter a{x), namely V ■ (ctVm) = 0, for the electrical potential u, the 
invisibility question has already been answered affirmatively in the case of 
electrostatics. 

The present work addresses the invisibility problem at all frequencies k ^ 0. 
We wish to cloak not just a passive object, but rather an active "device", 
interpreted as a collection of sources and sinks, or an internal current, within 
D. Since the boundary value problems in question may fail to have unique 
solutions (e.g., when —k"^ is a Dirichlet eigenvalue on D), it is natural, as in 
[GLUlj . to use the set of Cauchy data at dQ of all of the solutions, rather than 
the Dirichlet-to-Neumann operator on dQ, which may not be well-defined. 

The basic idea of |GLU3t [Lei IPSSlj is to form new EM material parameters 
by pushing forward old ones via a singular mapping F. Solutions of the rel- 
evant equations, Helmholtz or Maxwell, with respect to the old parameters 
then push forward to solutions of the modified equations with respect to the 
new parameters outside the cloaked region. However, when dealing with a 
singular mapping F, the converse is not in general true. This means that out- 
side D, depending upon the class of solutions considered, there are solutions 
to the equations with respect to the new parameters which are not the push 
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forwards of solutions to the equations with the old parameters. Furthermore, 
it is crucial that the solutions be dealt with on all of Q, and not only outside 
D. Especially when dealing with the cloaking of active devices, this gives rise 
to the question of what are the proper transmission conditions on dD, which 
allow arbitrary internal sources to be made invisible to an external observer. 
To address these issues rigorously, one needs to make a suitable choice of the 
class of weak solutions (on all of Q) to the singular equations being consid- 
ered. For both mathematical and, even more so, physical reasons, the weak 
solutions that are appropriate to consider seem to be the locally finite energy 
solutions; these belong to the Sobolev space with respect to the singular 
volume form \'g\^^'^dx on Vl for Helmholtz, and \'g\^^'^dx) for Maxwell. 

These considerations are absent from jLe| IPSSll IPSS2] , where the cloaking is 
justified by appealing to both the transformation of solutions on the exterior 
of D under smooth mappings F (essentially the chain rule) and the fact, 
in the high frequency limit, that rays originating in Q \ D avoid dD and 
do not enter D. As we will show, analysis of the transmission conditions 
at dD shows that the constructions of [PSSlt IPSS2j . although adequate for 
cloaking active devices in the absence of polarization, i.e., for Helmholtz, and 
cloaking passive devices in the presence of polarization, i.e., for Maxwell, fail 
to admit finite energy solutions to Maxwell when generic active devices are 
present. Furthermore, analysis of cloaking of an infinite cylinder, which was 
numerically explored in [CPSSPj and provides a model of the experimental 
verification of cloaking in [SMJCPSS] , shows that even cloaking a passive 
object may be problematic. Fortunately, it is possible to remedy the situation 
by augmenting or modifying this construction. 

We now describe the results of this paper. For what we call the single coat- 
ing, which is the construction of [GLU3j . and apparently that intended in 
|PSSlj . we establish invisibility with respect to the Helmholtz equation at 
all frequencies k ^ 0. In fact, one can not only cloak a passive object in 
a region D GG containing material with nonsingular index of refraction 
n{x), from all measurements made at the boundary dQ, but also an active 
device, interpreted as a collection of sources and sinks within D. 

Among the phenomena described here is that finite energy solutions to the 
single coating construction must satisfy certain "hidden" boundary condi- 
tions at dD. For the Helmholtz equation, this is the Neumann boundary 
condition at dD~ , and it follows that waves which propagate inside D and 
are incident to dD~ behave as if the boundary were perfectly reflecting. Thus, 
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the cloaking structure on the exterior of D produces a "virtual surface" at 
dD~. However, for Maxwell's equations with electric permittivity e{x) and 
magnetic permeability the situation is more complicated. The hidden 

boundary condition forces the tangential components of both the electric 
field E and magnetic field H to vanish on dD~ . For cloaking passive objects, 
for which the internal current J = 0, this condition can be satisfied, but for 
generic J, finite energy time-harmonic solutions fail to exist, and thus the 
single coating construction is insufficient for invisibility. In practice, even for 
cloaking passive objects, this may degrade the effective invisibility. 

We find two ways of dealing with this difficulty. One is to simply augment 
the single coating construction around a ball by adding a perfect electrical 
conductor (PEC) lining at dD, in order to make the object inside the coating 
material to appear like a passive object. (Such a lining was apparently incor- 
porated, although claimed to be unnecessary, into the code used in |CPSSPj 
in an effort to stabilize the numerics.) However, for the sake of brevity, the 
necessary weak formulation of the boundary value problem for this setup will 
not be considered in this paper. 

Alternatively, one can introduce a more elaborate construction, which we 
refer to as the double coating. Mathematically, this corresponds to a singular 
Riemannian metric which degenerates in the same way as one approaches dD 
from both sides; physically it would correspond to surrounding both the inner 
and outer surfaces of D with appropriately matched metamaterials. We show 
that, for the double coating, no lining is necessary and full invisibility holds 
for arbitrary active devices, at all nonzero frequencies, for both Helmholtz 
and Maxwell. It is even possible for the field to be identically zero outside 
of D while nonzero within D, and vice versa. 

Finally, we also analyze cloaking within an infinitely long cylinder, D C M'^. 
In the main result of §7] and ^ we show that the cylinder D becomes invisible 
at all frequencies if we use a double coating together with the so-called soft 
and hard (SHS) boundary condition on dD. For the origin and properties 
of the SHS condition and a description of how the SHS condition can be 
physically implemented, see |HLSl IKill IKi2t [Li] . 

We point out that there is some confusion in the physical literature concern- 
ing the theoretical possibility of invisibility. By this we mean uniqueness 
theorems for the inverse problem of recovering the electromagnetic param- 
eters from boundary information (near field) or scattering (far field) at a 
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single frequency, or for all frequencies. There is a vast literature on this 
subject. We only mention here mathematical results directly related to the 
one mentioned in [Let ISMJCPSS] . The Helmholtz operator at non-zero en- 
ergy for isotropic media is given by A -|- k^n{x), where n{x) is the index of 
refraction and /c 7^ 0. Unique determination of n{x) from boundary data for 
a single frequency k, under suitable regularity assumptions on n{x) and in 
dimension n > 3, was proved in |SyU] , with a similar result for the acoustic 
wave equation in [N]. (See [U] for a survey of related results). The article 
[N] was referred to in [Lei ISMJCPSS] as showing that perfect invisibility was 
not possible. However, the results of |SyU[ IN] for the Helmholtz equation 
are valid only under the assumption that the medium is isotropic and that 
the index of refraction is bounded. This does not contradict the possibility 
of invisibility for an anisotropic index of refraction, nor for an unbounded 
isotropic index of refraction. The constructions of jGLU3l iLel IPSS lJ and 
the present paper violate both of these conditions. We also point out that 
the counterexamples given in [GLUH Sec. 4] yield invisibility for the Helm- 
holtz equation, in dimension n > 3, for certain isotropic negative indices of 
refraction which are highly singular (and negative) on Q\D. 

We also note here that, at fixed energy, the Cauchy data is equivalent to 
the inverse scattering data. The connection between the fixed energy inverse 
scattering data, the Dirichlet-to-Neumann map and the Cauchy data is dis- 
cussed, for instance, in [N] for the Schrodinger equation and in [U] for the 
Helmholtz equation in anisotropic media. The scattering operator is well 
defined for the degenerate metrics defined here; see, e.g., [M]. 

There is a large literature (see [U]) on uniqueness in the Calderon problem 
for isotropic conductivities under the assumption of positive upper and lower 
bounds for a. It was noted by Luc Tartar (see [KVj for an account) that 
uniqueness fails badly if anisotropic tensors are allowed, since if F : Q — > Q 
is a smooth diffeomorphism with F\g^ = id, then F*cr and a have the same 
Dirichlet-to-Neumann map (and Cauchy data.) Note that since e and jj, 
transform in the same way, this already constitutes a form of invisibility, 
i.e., from the Cauchy data one cannot distinguish between the EM material 
parameter pairs e, /i and e = Jl = F^^jji. 

Thus, uniqueness for anisotropic media, in the mathematical literature, has 
come to mean uniqueness up to a pushforward by a (sufficiently regular) 
map F . Such uniqueness in the Calderon problem is known under various 
regularity assumptions on the anisotropic conductivity in two dimensions 
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[a EH [guUl lALP] and in three dimensions or higher [LiUl [LiUl ILTU] . 

but for all of these results it is assumed that the eigenvalues of a{x) are 
bounded below and above by positive constants. Related to the Calderon 
problem is the Gel'fand problem, which uses Cauchy data at all frequencies, 
rather than at a fixed one; for this problem, uniqueness results are also 
available, e.g., [BeKjlKK] , with a detailed exposition in [KKLj . For example, 
in the anisotropic inverse conductivity problem as above, Cauchy data at all 
frequencies determines the tensor up to a diffeomorphism F : Q — > Q. 

Thus, a key point in the current works on invisibility that allows one to 
avoid the known uniqueness theorems for the Calderon problem is the lack of 
positive lower and upper bounds on the eigenvalues of these symmetric tensor 
fields. In this paper, as in |GLU3i [Lei IPSSl] , the lower bound condition is 
violated near dD, and there fails to be a global diffeomorphism F relating 
the pairs of material parameters having the same Cauchy data. 

For Maxwell's equations, all of our constructions are made within the context 
of the permittivity and permeability tensors e and fi being conformal to 
each other, i.e., multiples of each other by a positive scalar function; this 
condition has been studied in detail in |KLS] . For Maxwell's equations in the 
time domain, this condition corresponds to polarization-independent wave 
velocity. In particular, all isotropic media are included in this category. This 
seemingly special condition arises naturally from our construction, since the 
pushforward {e, Jl) of an isotropic pair {e, fi) by a diffeomorphism need not 
be isotropic but does satisfy this conformality. For both mathematical and 
practical reasons, it would be very interesting to understand cloaking for 
general anisotropic materials in the absence of this assumption. 

We believe that our results suggest improvements which can be made in phys- 
ical implementations of cloaking. In the very recent experiment [SMJCPSS] , 
the configuration corresponds to a thin slice of of an infinite cylinder, inside 
of which a homogeneous, highly conducting disk was placed in order to be 
cloaked. This corresponds to the single coating with the metric g2 (see §2) 
on D being a constant multiple of the Euclidian metric. The analysis here 
suggests that lining the inside surface dD~ of the coating with a material 
implementing the SHS boundary condition |HLSl IKill IKi2l [Li] should result 
in less observable scattering than occurs without the SHS lining, improving 
the partial invisibility already observed. 

The paper is organized as follows. In §2 we describe the single and double 
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coating constructions. We then establish cloaking for the Helmholtz equation 
at all frequencies in §3. The notion of a finite energy solution for the single 
coating is defined in §§3.2 and then the key step for showing invisibility is 
Proposition 3.5. We discuss the Helmholtz equation for the double coating 
In §§3.3; there we define the notion of a weak solution and the Neumann 
boundary condition at the inner surface of the cloaked region. The key step 
for invisibility for Helmholtz at all frequencies in the presence of the double 
coating is Proposition 3.11. 

In §4 we study invisibility at all frequencies for Maxwell's equations. We 
define the notion of finite energy solutions for the single and double coat- 
ings. In §5 we demonstrate invisibility for Maxwell's at all frequencies for 
the double coating; see Proposition 5.1. In §6 we show that, for the single 
coating construction, the Cauchy data for Maxwell's equations must vanish 
on the surface of the cloaked region, showing that generically finite energy 
solutions for Maxwell's equations in the cloaked region do not exist. In §7 we 
consider an infinite cylindrical domain and show invisibility at all frequencies 
for Maxwell's equations for the double coating; the key result is Proposition 
7.1. In §8, we consider how to cloak the cylinder, treating its surface as an 
obstacle with the SHS boundary condition. Finally, in §9, we briefly indicate 
how general the constructions can be made. In particular, we note that a 
modification the double coating allows one to change the topology of the 
domain and yet maintain invisibility. 

We would like to thank Bob Kohn for bringing the papers [Lei IPSSlj to our 
attention, and Ismo Lindell for discussions concerning the SHS boundary 
condition. 

2 Geometry and basic constructions 

The material parameters of electromagnetism, namely the conductivity, (j{x)] 
electrical permittivity, e{x)] and magnetic permeability, fi{x), all transform 
as a product of a contravariant symmetric 2-tensor and a (+1)— density. That 
is, if F : Qi — > Q2, y = F{x), is a diffeomorphism between domains in 
M", then a{x) = (cr-'^(x)) on Qi pushes forward to {F^a){y) on Q2, given by 
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1 




(1) 



p,q=l 



x=F-^(y) 



with the same transformation rule for the other material parameters. It was 
observed by Luc Tartar (se |KV] ) that it follows that if F is a diffeomorphism 
of a domain fixing dVL, then a and a := F^,a have the same Dirichlet-to- 
Neumann map, producing infinite-dimensional families of indistinguishable 
conductivities. 

On the other hand, a Riemannian metric g = {gjk{x)) is a covariant sym- 
metric two-tensor. Remarkably, in dimension three or higher, a material pa- 
rameter tensor and a Riemannian metric can be associated with each other 



where {g^^) = {gjk)~^ and \g\ = det {g). Using this correspondence, examples 
of singular anisotropic conductivities in M",n > 3, that are indistinguishable 
from a constant isotropic conductivity, in that they have the same Dirichlet- 
to- Neumann map, were given in |GLU3] . The two constructions there are 
based on two different types degenerations of Riemannian metrics, whose 
singular limits can be considered as coming from singular changes of vari- 
ables. The singular conductivities arising from these metrics via the above 
correspondence are then indistinguishable from a constant isotropic a. In 
the current paper, we will continue to examine one of these constructions, 
correpsonnding to pinching off a neck of a Riemannian manifold; we refer to 
it as the single coating. We also introduce another construction, referred to 
as the double coating. We start by giving basic examples of each of these. 

For both examples, let Q = -8(0,2) C M^, the ball of radius 2 and center 
0, be the domain at the boundary of which we make our observations; D = 
-5(0, 1) C the region to be cloaked; and S = dD = EP' the boundary of the 
cloaked region. 

Single coating construction: We begin by recalling an example from 
|GLU3l [PSST] ; the two dimensional examples in [Le| IV] are either essentially 
the same or closely related in structure. 



by 




or g^'' = \a 



2/{n-2)ik 



(2) 
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For the single coating, we blow up using the map 

Fi : S(0,2)\{0} A = (^ + 1)-, r = < r < 2. (3) 

On B{0, 2), let {ge)ij = ^ij be the Euclidian metric, corresponding to constant 
isotropic material parameters; via the map Fi, pushes forward, i.e., pulls 
back by F^^, to a metric on Q\D, 

g, = (Fi).^?e := (F-'Yig,) . 

Introducing the boundary normal coordinates {u, r) in the annulus Q\ D, 
where u = {u^, u"^) are local coordinates on S = and r > is the distance 
in metric gi to S, we have, from ([3]), 

gi = du^ + dr'^, r = 2(r - 1). (4) 

Here duj"^ = hap{uj)duj'^duj^ is the standard metric on induced by the 
Euclidian metric on M?. Note that 'gi has the following properties: 

Consider a local (yfe-orthonormal frame {driV^w) oiiVL\D consisting of the 
radial vector 

^ d d 
^ dr r dx^ 
and two vector fields v,w. Then, 

gi{dr,dr) = gi{dr,v) = gi{dr,w) = 0, gi{w,v) = 0, (5) 

_ 1)2 e [Cl' C2], e [Ci, C2], 

where ci,C2 > 0. Thus, 'gi has one eigenvalue bounded from below (with 
eigenvector corresponding to the radial direction) and two eigenvalues that 
are of order (r — 1)^ (with eigenspace span{v^ w}). In Euclidean coordinates, 
we have that, for \gi\ = det (^^i), 

-C,{r-l)\ (6) 
2x 

Wi^il <C2, T^i = — = 2(<9r)i. 

r 

Here and below we use Einstein's summation convention, summing over 
indices appearing both as sub- and super-indices in formulae, and u = 
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(z/i, z/2, z^s) denotes the unit co- normal vectors of the surfaces {x & D,\ D : 
|a;| = s}, l<s<2, with respect to the metric g. 

On D, we simply let 'g2 be the Euclidian metric. Together, the pair (^1,^2) 
define a singular Riemannian metric on fl, 



which is singular on E+, i.e., as one approaches E from D,\D; in the sequel, 
we will identify the metric g and the corresponding pair {gi,g2)- 

To unify notation for the two basic constructions, we will denote in the single 
coating case Mi = Q, M2 = D and let M be the disjoint union M = M1UM2. 
Also, for notational unity with the double coating, we let 71 = {0} C Mi, 
72 = C M2, and 7 = 71 U 72. Moreover, we denote Ni — il \ D, N2 = D, 
E = and = A^i U E U A^2 := C R^. 

Double coating construction: The double coating refers to a metric on Q 
that is degenerate on both sides of E and has the same limit as one approaches 
S from both sides. 

We now introduce notation, shared with the single coating, that will be used 
throughout for the double coating. Let Mi = Q = B{0,2), which is compact 
with boundary, and M2 := S'^^^, the 3-sphere of radius I/tt, which is compact 
without boundary, and again let M — M1UM2 be their disjoint union. For the 
double coating, let 71 = {0} C Mi, 72 = {NP} C M2, where NP is a chosen 
point, e.g., the North Pole of S'l/^, and 7 = 71 U 72. As in the single coating 
example, we let A^i = n\D = B{0,2)\B{0, 1), N2 = D = B{0, 1), S = dD, 
and A^ = A^i U S U A^2 C M^. We take the diffeomorphism Fi : Mi \7i — > Ni 
to be the blow-up of 71 as in the single coating, while we blow-up 72 by defin- 
ing F2 : M2\72 — A^2 as follows. Denote by SP the point on Q2 antipodal to 
NP. Then the Riemannian normal coordinates centered at SP are defined on 



Introduce (local) spherical coordinates {00, r) on A^2 = B{0, 1), considered as 
a subset of Tsp{E>^), with u; = (u;i,a;2), G E = dB{0, 1) and < r < 1. 




xen\D, 

X e D, 



B{0, 1) C TspE>^ ~ 

exp^p: B(0,1)^M2\{NP}. 
Denote by F2 the diffeomorphism 



F2 = (exp5p)-^ : M2 \ {NP} ^ B{0, 1). 
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The standard metric g on S^^^ in these coordinates takes the form 

sin^fvrr) ,072 /w\ 
92 = duj^ + dr\ (7) 

where duo"^ is again the standard metric on 

Observe that g2 = (-p2)*(fi'2), as one approaches S~ on 5(0,1), has very 
similar properties to gi on 5(0,2) \ B{0, 1) as one approaches S+. Indeed, 
again consider the radial vector dr = ^ = ^gfj aX x E N2 and two vectors 
f,w such that in Euclidean metric {dr,v,w) is a local orthonormal frame. 
Then, as follows from (^^, ai x & N2 with, say, 1/2 < r < 1, 



92idr, dr) = 1, 5'2(<9r, ^^) = fi'2(c^r, ti?) = 

:i-r)2' fl 



^2(w, w) = 0, _ G [ci, C2] 



where ci,C2 > 0. Thus, g2 has one eigenvalue bounded from below (with 
eigenvector corresponding to the radial direction) and two eigenvalues that 
are of order (1 — r)^, and with respect to the Euclidean coordinates on N2, 

\92{x)\'^' < Ci(l - r)^ \gl'u.\ < C2, v, = -- = -(9,)„ ^ < r < 1. (8) 

r 2 

Set 'gi = {Fi)^ge on A''i, where Fi is defined as for the single coating example. 
Together, these define a singular metric g = (gi,g2) on the entire ball N = 
A^i U A^2 U S = 5(0,2). Comparing (HD and ([7]), we see that, in the Fermi 
coordinates associated to S, \g\^^'^'g'^-' is Lipschitz continuous on A^; note 
also that \'g\^^'^])ij is not invertible at dB{0, 1). 

Although they are distinct, each of these constructions may be summarized 
as follows. The domain Q, which we will refer to as A^, decomposes as 
A^ = A^i U S U A^2, where Ni = \D, N2 = D and J: = dD. Ni and 
A'2 are manifolds with boundary, with dNi = dQ U dD~^ = dN U S+ and 
dN2 = E~, where the superscripts ± are used when considering limits from 
the exterior or interior of the cloaked region. The singular electromagnetic 
material parameters on A^ will correspond to a singular Riemannian metric 



^Recall that the Fermi coordinates associated to S are {uj,t), where ui — {uj^,uj'^) are 
local coordinates on E and r = t{x) is the distance from a; to S with respect to the metric 
multiplied by +1 in A^i and —1 in N2- 
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g = {gi, g2), arising as the pushforward of a (nonsingular) Riemannian metric 
g = (fl'i)5'2) on a manifold with two components, M = Mi U M2, via a map 
F -.MX-i — ^ N, 

^ ' \ F^ix), a; G M2\7. 

Here, Mi and M2 are disjoint, with Mi diffeomorphic to iV; 71 = 7 fl Mi is 
either a point (the point being blown up) for the single and double coatings, 
or a line (for the cloaking of an infinite cylinder in ^ }7f8|) ; and 72 = 7 H M2 is 
either empty (for the single coating) or a point (for the double coating) or a 
line (for the cylinder.) In §21 we will show that such constructions exist in 
great generality, and for this reason the proofs will be expressed in terms of 
analysis on M and A^. 

In this generality, we say that [M, N, F,'-f,J], g) is a coating construction if 
(M, g) is a (nonsingular) Riemannian manifold, 7 C M and S C are as 
above, and F: M\'y^N\'Eis diffeomorphism of either type. This then 
defines a singular Riemannian metric g everywhere on \ E = A^i U A''2 , by 

f £1 := Fugi, X G A^i, 
\g2-= ^2*^2, X G N2. 

If we introduce Fermi coordinates {u, r) near S as above, the g satisfies 
([S]),([HD or ([H]), with r — 1 replaced by r, for the single and double coatings, 
resp. From these, one sees that \g\^^'^g-''' has a jump discontinuity across S 
for the single coating and is Lipschitz for the double coating. Note that in 
both examples, N = Q = B{0,2), so that and Mi have the same topology. 
However, in a direct extension of the double coating construction, described 
in §ni the domain A^ containing the cloaked region A'2 need not even be 
diffeomorphic to Mi. 



3 The Helmholtz equation 

We are interested in invisibility of a cloaked region with respect to the Cauchy 
data of solutions of the Helmholtz equation, 

(Ag + k'^)u = f in (9) 
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where / represents a collection of sources and sinks. The Cauchy data ^ 
consists of the set of pairs of boundary measurements {u\qq,, dyu\Qa) where u 
ranges over solutions to ([9]) in some function or distribution space (discussed 
below). Let (M, A^, F, 7, S, g) be a single coating construction as in ^ For 
the moment, as in the Introduction, we continue to refer to as fi, N2 as 
D and S+ as dD^\ we may assume that Mi = N, M2 = D and F2 = id, 
so that g2 = (72 is a (nonsingular) Riemannian metric on D. Thus, 5^ is a 
Riemannian metric on Q, singular on Q\ D, resulting from blowing up the 
metric gi on Q with respect to a point O and inserting the {D,g2) into the 
resulting "hole". 

We wish to show that C| j = Cg,o frequencies < < 00, if supp(/) C 

D and k is not a Neumann eigenvalue of {D,g2). Due to the singularity of 
g, it is necessary to consider nonclassical solutions to (Q, and we will see 
that the exact notion of weak solution is crucial. Furthermore, a hidden 
Neumann boundary condition on dD~ is required for the existence of finite 
energy solutions. Physically, this means that the coating onQ\D makes the 
inner boundary dD~ appear to be a perfectly reflecting "sound-hard surface" 
for waves propagating in D, while, from the exterior, the cloaked device is 
invisible; that is, measurements of solutions of the Helmholtz equation at dQ 
cannot distinguish between {Q,g) and (^,(7). 

3.1 k = and weak solutions 

First consider the case when k = and / = 0. As described in the In- 
troduction, this situation was treated in [GLU3j in the context of electrical 
impedance tomography. There, it sufficed to consider as weak solutions those 
L°° functions satisfying (for the metric g) in the sense of distributions. It 
was shown that, for given Dirichlet data h on dfl, ([9]) has a unique such solu- 
tion, u, which must, by removable singularity considerations, be constant on 
D. These same conclusions would have held if we had considered the larger 
class of spatial weak solutions (defined below). However, for > or 
/ 7^ 0, we will see that this notion of weak solution is inappropriate. 
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3.1.1 A; > and spatial solutions 

Definition 3.1 u is a spatial solution to the Dirichlet problem for the 
Helmholtz equation, 

{Ag + k'^)u = f onQ,u\9n = h (10) 

^f 

ueH\n,dx) (11) 

and 

dM'^'rd,u) + e\g\'/'u=\g\'/'f m H-\^,dx). (12) 

Here, for s G M, H^{Q,dx) refers to the standard Sobolev space of distribu- 
tions with s derivatives in L'^{Q, dx). Note that flTT]) . together with the prop- 
erties of the metric tensor given in ^ implies that \g\^^'^g^^diU G L^(i7, dx). 

Later in our analysis (see 0361) ). we will see that fll2p implies that the normal 
derivative of u on dD~ vanishes, 

dru\dD- = 0. 

On the other hand, the fact that u G H^{Q, dx) implies that 

u\dD- = u\qd+ = constant := u{0), 

with u the solution to (Ag + k'^)u = in dQ, uIqq = h, where the first 
equality follows from the trace theorem for functions and the second 
from considerations similar to those in [GLUSt Prop. 1]. Note that, for 
generic k and h, u{0) ^ 0. Thus, U2 := m|d needs to be a solution of the 
overdetermined elliptic boundary value problem on {D,g2), 

(A + k'^)u2 = 0, di,U2\dD = 0, U2\dD = constant 7^ 0. (13) 

Clearly, for generic A; > there exists no solution to f|T3|) and therefore there 
is no weak solution to ffTOl) in the sense of Definition 13. II Rather, one needs to 
use an norm adapted to the singular Riemannian metric g; this is in fact 
physically natural, being essentially the energy of the wave. We formulate 
the correct notion in the next section. 
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3.2 Finite energy solutions for the single coating 



We now give a more satisfactory definition of weak solution, restricting the 
notion to those solutions that are physically meaningful in that they have 
finite energy. 

We now revert to the notation of M, N, . . . when discussing the single coat- 
ing construction, i.e., let {M, N, F,'y,J2, g) denote a single coating as in ^ 
Our first task is to understand in what sense the expression \'g\^^'^'g^^diU is 
rigorously defined. 

To this end, define for G C°°(N) 

J N 

Let 

be the completion of C°°{N) with respect to the norm || ■ We note that 
H^(N, \g\^^'^dx) C L'^{N, \'g\^^'^dx), so we can consider its elements as mea- 
surable functions on N. 

Lemma 3.2 The map 

has a bounded extension 

Dg : H^{N, \g\^''^dx) M{N; M^), 

where Ai[N;M.^) denotes the space of M^^-valued signed Borel measures on 
N. Moreover, for u ^ X , we have, in the sense of Borel measures 

(D^«)(S) = 0. (14) 

Proof. Let G C~(iV) and ^ G C(N). Then Dl^) G L°°{N). Let = 
(f)oF,ri = rjoFe Then, 

/ {D~(j))rjdx = / {D~(j))rjdx 
Jn Jn\t, 

= [ \9\'^'9''^d,c^vdx+ [ \g\'/V'%d,c^vdx. 
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As the metric g is bounded from above and below, and ^ = 0{r ^) on Mi 
and = 6j on M2, we have 



{D-d)) 7] dx 



N 



< Co{\\(t)\\m(Mx,dx)\\'n/A\L-^{Mi4x) + \\<P\\m{M2,dx)\\v\\L-2{M2 4x)) 



< C'i||0lUI|r/||c(iv)4^^(supp(?7),S) 



where dg is the distance on N with respect to the metric g. This shows the 
existence of the bounded extension Dg : H-^{N, \g\-^^'^dx) — > M{N; R^). Also, 
if we consider functions rj supported in small neighborhoods of S, we see that 
([HD follows. □ 

We also need the following auxiliary result 

Lemma 3.3 Assume that u is a measurable function on N such that 

u e L^{N, \g\^''^dx), (15) 
u\N\i:&Hlc{N\T.,dx), (16) 

/ \g\^^^g'^diudjudx < 00. (17) 



Then u G H^{N, \g\^/'^dx). 



Note that, due to the fact that 'g is bounded and positive definite on any 
compact subset of \ S, condition f|T6|) in fact follows from conditions ([T5|) . 
(fTTl) and is included for the convenience of future references. 



Proof. Consider first the case when m = in A^i. 

First, the condition f|T7j) implies that v = u\n2 ^ H^{N2,dx). Let / = f |s G 
H^/'^iE) and G H\Ni, dx) be an extension of /. Let x e C^(K) be a cut- 
off function with x(t) = 1 for \t\ < | and x(^) = for \t\ > 1. We introduce 
Fermi coordinates near S as in §2], (r, tu), r G (0,2), u = (a;i,ti;2) G S. 

Define, for e > 0, 

/ \ _ / v{x), X G A'2, 
"^^^"^^ \ X(^)i5;^(a:), x G ATi. 

Then G H^{N,dx) and, using ([3]), ([5]), we see that 

lim / \g\^/'^[g'^di{we-u)dj{we-u) + {we-uy]dx (18) 

= lim / + = 0, 

Jni 
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Observe that the integrand vanishes outside the a neighborhood of S"*" of 
volume less than Ce. Next, divide the integral involving derivatives in the 
right-hand side of f|T8|) into the terms involving components tangential and 
normal to the boundary, using the fact that r = 2(r — 1): 

JNi\T, ^ 

and where a, /3 run over {1, 2}, 



1^ 



1/2 



dr[xC-)Ef] 



e 



dTduJiduj2- 



As, by ([5]), \'g\'^^'^'g"'^ is bounded, the integral involving tangential derivatives 
tends to due to the volume of the domain of integration. Again, by ([5]) 
we have {g]^^"^ ^ Ct"^; this, together with the volume estimate and the fact 
that |c}^x(~)| ^ Ct~^, implies that the integral involving normal derivatives 
tends to when £ — > 0. Similarly, we see that 



/ \g\^^^\xi-)E^\^dx -^0 forego. 



The function G {N, dx) can be approximated with an arbitrarily small 
error in H^{N,dx) by a C°°{N) function, and we see that the same holds 
in the X-norm. Thus We G H^(N, \'g\^/'^dx), and the above limit shows that 

ueH\N, \g\^/^dx). 

Now let M be a measurable function in satisfying f|T5|) . f|T6l) . and f|T7|) . 
Let be the characteristic function of A^2- As XN2U G H^iN, \g\^^'^dx), it 
is enough to show that u — XN2U G H^{N, \g\^/'^dx). This means that it is 
enough to consider the case when m = in A^2- Clearly, we can restrict our 
attention to the case when u vanishes also near dN. 

Now let ui = u o F in Mi \ 71. Then we see that 



\'/^g'^d,{u,)d,{u,)dx< 



00. 



MiVa 



Let w = VuImiX-ji- Using a change of coordinates in integration and f[T5|) . 
we see that u G L^(Mi \ 'yi,dx). Extending ui and w to functions ul and 
on 7i, we obtain functions ul G L'^{Mi,dx) and M^-valued function 
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G L'^{Mi,dx). Now Vm^ — E H^^{Mi, dx) is supported on 71. Since 
there are no non-zero H~^{Mi, dx) distributions supported on 71, we see that 
'Wu\ = w'^ E L^(Mi, dx). Thus we see that G H^{Mi, dx). In the following 
we identify ui and ul. As ui vanishes near dMi, and 71 consists of a single 
point and thus is a (2, l)-polar set [Maj pp. 393-7], there are (pj G {Mi\ji) 
such that (pj —>■ Ui in H^{Mi, dx) as j 00, that is, 

lim / \g\^%'''di{^j - u)dki(f)j - u) + - u)^] dx = 0. 

J Ah 

Now let 4>j G C^{N), with supp(0j) C A^i and 

~ _ r 0,oFri iniVi, 
\ iniVa. 

Then the previous equation implies that 

lim [ Igl'^^ig^'dSj - u)dSj - + i^j ' ^)^] dx = 
lim / Igl'^^ig^'dii^j - u)dk{4>j -u) + (0,- - uf] dx = 0, 

where we use that u = in N2. 

This shows that (pj is a sequence converging in the X-norm and that the 
limit is u. Thus u G H^(N, \'g\^^'^dx), proving the claim. □. 

Although in this section (M, N, F, 7, E, g) continues to denote a single coat- 
ing, we will see later that the following definition is also appropriate for the 
double coating construction. 

Let / G L'^{N, dx) be a function such that supp (/) fl S = 0. 

Definition 3.4 Let {M, N, F,'y,J], g) be a coating construction. A measur- 
able function u on N is a finite energy solution of the Dirichlet problem for 
the Helmholtz equation on N, 

{Ag + k'^)u = f onN, (19) 
u\dN = h, 
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ueL'^iN, \g\^^^dx); (20) 
u\N\i:eHUN\J:,dxy, (21) 

f \g\^/Y^diudjudx < oo, (22) 
Jn\j: 

u\dN = h, ; 
and, for all ^ G C°°{N) with ^{g^ = 0, 

-{Dgu)dji) + k'^uij\g\^''^]dx = [ f{x)^{x)\g\^^^dx (23) 

N Jn 

where the integral on the left hand side of [2^) is defined by distribution-test 
function duality. 

Note as before that condition follows from fl20|l . (|22D. 

Cloaking by the single coating of arbitrary active devices, with respect to 
solutions of the Helmholtz equation at all frequencies, then follows from the 
following. 



Theorem 3.5 Let u = (^1,^2) : M \ ^ R and u : N \ T. ^ R be 

measurable functions such that u = u o F . Let f = (/i, /2) : M \ 7 — M 
and f : N \ Ti —>■ R be L"^ functions supported away from 7 and S such that 
f = foF, andh: ON ^R, h: dMi -^Rbe such that h = hoFi. 

Then the following are equivalent: 

1. The function u, considered as a measurable function on N , is a finite 
energy solution to the Helmholtz equation / fTPj) with inhomogeneity f 
and Dirichlet data h in the sense of Definition \3.4 ■ 

2. The function u satisfies 

(Ag + k'^)ui = /i on Ml, MilaMi = K (24) 

and 

iAg + e)u2 = f2 onM2, g^'iyAu2\dM,=b, (25) 

with 6 = 0. Here ui denotes the continuous extension of ui from Mi \7 
to Ml 
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Moreover, ifu solves [24\ ) and l[2^) with b ^ 0, then the function u = uoF~^ : 
N\'E ^ M, considered as a measurable function on N, is not a finite energy 
solution to the Helmholtz equation. 

Remarks, (i) It follows that the construction of |GLUlt [PSSlj cloaks active 
devices from detection by unpolarized EM waves at all frequencies. 

(ii) Observe that in fl2^ the right hand side /i is zero near 71. Thus ui, 
considered as a distribution in a neighborhood of 71, has an extension on 71 
that is smooth function in a neighborhood of 71. 

(iii) As noted previously, for the single coating case one may assume that 
A^2 = M2 and F\m2 is the identity. Thus uIn^ = u\m2] hence, if m is a finite 
energy solution of the Helmholtz equation on A^, we see that u\m2 satisfies the 
Neumann boundary condition on dM2 and thus also 'u|7V2 automatically has 
to satisfy the Neumann condition on S~. The Neumann boundary condition 
that appears on dN2 means that, observed from the inside of the cloaked 
region N2, the single coating construction has the effect of creating a virtual 
sound hard, i.e., perfectly reflecting, surface at S. Similarly, we will see later 
that there are hidden boundary conditions for Maxwell's equations in the 
presence of the single coating, but they are overdetermined and generally 
preclude such solutions existing. 

Proof. First we proof that Helmholtz on M implies Helmholtz on A^. 

Let / e L'^{M, dx) be a function such that supp (/) n (7 U dMi U SMa) = 0. 
Assume that a function m on M is a classical solution of (12^ and fl25l) . Notice 
that we have required here that U2 on dM2 satisfies the Neumann boundary 
condition at dM2. 

Again, define u = F^u and f = f o F^^ on N \ T, and extend it, e.g., by 
setting it equal to zero on S. Note that then / G L'^{N,dx) is supported 
away from S, and u G L'^{N, \'g\^^'^dx) satisfies 

{Ag + k'')ui = fi = f\N, mNi, u\dN = h, (26) 

and 

{Ag + e)u2 = h = f\N2 inA^a. (27) 

Let ^(e) be the e-neighborhood of S with respect to the metric g. Let 7(e) 
be the e-neighborhood of 7 C Mi with respect to the metric g. Let gbnd and 
gbnd be the induced metrics on d'y{e) and dT,{e), correspondingly. 
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Clearly, the function u satisfies conditions (12U]) . and (12^ . By Lemma 
13.31 we have that u E H^{N, \g\^^'^dx), and DgU is thus well defined. 

Using relations for the normal component and fl2^ . flTTI) . and property 
([Till of DgU, we see that, for ^ G Co~(A^), 

/ [-/^^(n)^,^^ + eui^\g\'/^ - fm'^']dx (28) 
Jn 

= lim / {-g'^ diudji) + {k'^u-J)^)\g\^''^dx 
= \im{[ +[ ){-g'^u,d,ui^)\gbnd\'^^dS 

JdT.(e)nN2 JdT,{e)r\Ni 

= lim/ (-5^*Vi9,M2t^)|5^6„d|^/'t/S+ (29) 

+ lim / {-g^^ d,u, o F))\gbnX^'dS (30) 
= 0. 

Indeed, the integral (l29l) in the right-hand side of this equation tends to 
due to the boundary condition on ( [25|) . and boundedness of o F. To 
analyze the integral (!30!) observe that, as supp/i fl 71 = 0, mi is infinitely 
smooth near 71. Thus all diUi and ipo F are bounded near 71, while the area 
of 97(e) is bounded by Ce"^. Hence we see that (!23l) is valid and thus 

{Ag + k'^)u = f inN 

in the sense of the Definition 13. 4[ 

Summarizing, so far we have proven that a (classical) solution to the Helm- 
holtz equation on M yields, via the pushforward, a finite energy solution to 
the equation on A^. 

Next we consider the other direction and prove that the Helmholtz equation 
on N implies Helmholtz equation on M. 

Assume that u satisfies Helmholtz equation f|T9|) on {N,g) in the sense of 
Definition 13. 4[ with / G L?'{N) supported away from S. In particular, u is a 
measurable function in satisfying f|T5|) . f|T6l) . and ffTTj) . 

Let u = uo F and f = foF on M\'-f. Then we have 

(Ag + k'^)ui = fi = /U/i\7i in Ml \ 71, Mi|aMi = h (31) 
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and 

{Ag + e)u2 = f2 = flu, inMs. (32) 
By conditions (fT5|) . (fT6l) . and f|T7|) . we have that 

g^\d,u){dku) G Li(Mi \7i, l^^l^/'rfx). 

and thus ui G H^{Mi \ 71, cix). As before, we see that 

{Ag + P)ui = /i in Ml, MilaAfi = h, (33) 

where /i is extends to have the value at 71 and Ui is smooth near 71. 
Let us now consider the boundary condition on M2. Since u satisfies 
we see that for ijj e C^{N), 



= [-Dgud,ij + euij\g\'/^~fij\g)\'/^]dx (34) 
Jn 

= Urn {-f^ diudj^+{k'^u-f)i))\g\^'^dx 



9E{e)nAf2 JaT.{e)r\Ni 



hm( / + / ){-g'^ V, d,u^)\%^d\^l''dS 



\ {-9'^ '^jdiU2\dM2^)\9bnd\^^'^dS (35) 

JdM2 

+ hm / {-g'^^d,u,u^^)\gU'^'dS 

JdMe) 



'-g'^ Vjdiu\QM2'4')\9bnd\^^'^dS, 



where ip = ipoF. Here we use the fact that Ui is a smooth function, implying 
that diUi is bounded and that ip = ip o F is bounded. As ■iIj\qm2 ^ C°°{dM2) 
is arbitrary, this shows that 

g'^ yjdiU2\dM2 = 0. (36) 

Thus, we have shown that the function m is a classical solution on M of 

(Ag + fc2)ui = /i in Ml, u^\9M, = h (37) 
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and 

{Ag + e)u2 = f2 inM2, g'''ujdkU2\aM, = 0. (38) 
This proves the claim, and finishes the proof of Theorem I3.5[ □ 

3.2.1 Operator theoretic definition of the Helmholtz equation 

It is standard in quantum physics that a self-adjoint operator can be defined 
via the quadratic form corresponding to energy. In the case considered here, 
the energy associated with the wave operator is defined by the quadratic 
(Dirichlet) form A, 

A[u,u\:= [ \g\^/^g'^diudjudx, ueV{A) (39) 

As we deal with the sound-soft boundary dN or, more generally, with the 
source on c^A^ of the form u\sn = h, the domain T'(A) of the form A should 
be taken as 

V{A) = Hl^{N, \g\^^^dx) C X. 

Thus, by standard techniques of operator theory, e.g., |Kaj . the form A de- 
fines a positive selfadjoint operator, denoted Aq = — A~ , on L'^{N, {gl^^^dx). 
Next we recall this construction. We say that u G Hq(N, \'g\^^'^dx) is in the 
domain of Aq, u G ^{Aq) if there is an / G L?'{N, \g\^^'^dx) such that for all 
veH^iN, \g\^/^dx), 

A[u,lJ\ = [ fv\g\^/^dx. (40) 

In this case, we define 

AqU = J. 

Proposition 3.6 Assume that — /c^ is not in the spectrum of A~ . Then u 
is a finite energy solution to 

{Ag + P)m = 7, u\qj^ = he H'/\dN) 
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if and only if 

u = Eh + {Af + k^r\f- {Ag + e)Eh), (41) 
where Eh is an {N , dx) - extension of h to N satisfying supp{Eh) C dN U 

Proof. First we show that if u satisfies the conditions of Definition 13.41 then 
it satisfies dH]). As -0 e C^{N), Vlajv = 0, imply that ijj e H^iN, \g\^/^dx), 
we see by fl23l) that u — Eh satisfies 

[ {-Di{u~ Eh)djV + k^{u- Eh)^dx = I \g\^/\f-{Ag + k^)Eh)vdx, 

for any v E C^{N). By ([H]) and ([22]), this imphes 

/ Is-r^^ (-g'Wi{u - Eh)djV + k^{u- Eh)v) dx (42) 

= I \g\''\J-{Ag + e)Eh)vdx, 
Jn\t, 

for any v G C^{N). We need to show that fj42l) is vahd for all v e 
H^{N,\g\^/'dx). 

Observe that 

/ Ifl^r^^ (-g'^di{Eh)djV + k\Eh)v) dx= [ \g\^^\{Ag + e)Eh)vdx, 

iiV\E ^ ^ Jn\T. 

where we use that supp{Eh) C dN U Ni and v\dN = 0. Thus, it remains to 
show that 

[ \g\^'^ {-g'^diudjv + k'^uv) dx= I \g\^'^Jvdx (43) 
Jn\t. Jn\t. 

for V G Hq(N, \g\^^'^dx). Clearly, to show this it is enough to show that 

lim / \g\^^^ (-g'^diudjV + k^uv - fv) dx = 0. (44) 
where Y.^{e) = A^i n S(e). 
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Next we argue analogously to the reasoning that led to equation 0281) . Let 

v = voF,f = foF, and -u = m o F in M \ 7. To clarify notations, denote 

ui = m|mi, U2 = u\m2 , vi = v\mj_, V2 = v\m2, and /i = /|mi, /2 = /Ufa- 
Then, by Proposition I3.5[ 

(A3 + PK = /i, in Ml, (45) 

{Ag + e)u2 = /2, in M2, (46) 

duU2\dM2 = 0, (47) 

and we see that 



^^°i(Mi\7(e))UAf2 

g'^uj d,{u) v)\gU'^''dS + [ {-g'^u, d,iu) v)\gU'^^dS. 

J a Mo 



By (jUD, we have that 



f {-g'^u, d,{u) v)\g,nd\^'HS = 0. (48) 
JaM2 



Next we consider 



h{s)= / {-g'^u,d,{u)v)\gU^'^dS. 

Jd'y{e)nAh 



Note that linie^o h i^) exists as the limits fl44l) and (1481) exists. 

As supp (/) n 7 = 0, we see that ui is smooth function near 7. Moreover, as 
■u G X, we observe that Vi G H^{Mi\'j, dx), and so vi can be extended to vi G 
H^{Mi,dx). Hence, by the Sobolev embedding theorem, vi G L^{Mi,dx). 
This allows us to deduce that 



liminfe-3/2 / \vi\ dS = 0. (49) 



Indeed, 



\vi\dS{x) \ dr = / |t>i|(ia; 
<( [ Ivil'^dx] ' ( f da;V =o(£^/'). 



7(e) / \Jl{e) 
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Clearly, this inequality implies (H^ . Thus using boundedness of g^^ujdiU we 
see that 

liminf [ {-g'hj di{u) v\gbnd\^'^)dS = 0. 

As \\me^Qli{e) exists, this implies \\m^_,Qli{e) = 0. As u\g]y = hhj Defini- 
tion [33] we have shown that Definition 13.41 implies fj4T|) . 

Next, consider the case when u satisfies (HTil . Since m G X, we see by ( |T^ 
that 

/ Di{u)djvdx= [ \g\^/^g'^diudjvdx (50) 

for all V e <:7o°°(A^). Thus, by gl]) we see that (gSD is valid for v G Co°°(A^), 
which implies condition (I22p . The other conditions in Definition 13.41 follow 
easily from fl41l) . □ 



3.3 Helmholtz for the double coating 

We now examine solutions to the Helmholtz equation in the presence of the 
double coating; we will establish full-wave invisibility at all nonzero frequen- 
cies. Unlike for the single coating, for the double coating no extra boundary 
conditions appear at S. Otherwise, the reasoning here parallels that in §3.2. 

Throughout this section, {M, N, F,'y,J2, g) is a double coating construction. 
3.3.1 Weak solutions for the double coating. 

Suppose that A; > and / G L'^{N, \g\^^^dx). We use the same notion of weak 
solution as for the single coating, saying that -u is a finite energy solution of 

{Ag + k^)u = f inN, u\qj^ = h (51) 

if u is a solution of the Dirichlet problem in the sense of Definition 13. 4[ 

We start with analogues of the space H^{N, \'g\^^'^dx), and Lemmas 13.21 and 
1331 To this end define, for G C°°{N), 

Il0lly := / i\g\'^'rd,^d,^+\g\'/'mdx. 
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Let 

H\N, \g\^/^dx) =Y := cly(C°°(iV)) 

be the completion of C°°{N) with respect to the norm || ■ ||y. Note that 
H^{N, I'gl^^'^dx) C L'^{N, \'g\^^'^dx), so we can consider its elements as mea- 
surable functions in N. 

Lemma 3.7 The map 

has a bounded extension 

Dg : H\N, \g\^^^dx) M{N; M.^), 

where M.{N] M'^) denotes the space of M^^ -valued signed Borel measures on 
N . Moreover, for u , we have 

{Dgu){i:) = 0. (52) 

Ifuis a measurable function on N such that 

u e L\N, \g\^/^dx), (53) 
M|iv\E e Hl^{N \ S, dx), and (54) 

/ \g\^^'^g^^ diudjU dx < oo, (55) 
Jn\t, 

then u e H^{N, \g\^/^dx). 

Proof. The proof here is essentially the same as of Lemmas 13.21 and 13.31 
The only difference is that, as described in ^ the map 

now consists of two maps, 

F,:MA7^iV„ 2 = 1,2, 

having similar structure to each other, namely that of the map Fi in the 
single coating construction. (Recall that for the double coating construction, 
7i := 7 n Ml is a point O E Mi and 72 := 7 n Ma a point NP E Mg.) 



29 



Therefore, when proving that u satisfying fl53l) - fl55l) is in H^(N, \'g\^^'^dx), we 
can use the fact that, in this case, both (1 — Xa^i)^ and (1 — Xn2)u satisfy 
f l53|) - fl55|) and carry out the proof for each of them as for the (1 — Xn2)'^ term 
in the proof of Lemma I3.3[ 

Invisibihty of active devices in the presence of the double coating with respect 
to the Helmholtz equation at all frequencies then follows from 

Theorem 3.8 Letu = (^1,^2) : M\7 M andu : A^\S ^ M 6e measurable 
functions such that u = uoF . Let f = (/i, : M\7 M and f : A^\S —>■ M 
be LP' functions supported away from 7 and S such that f = f o F. Then the 
following are equivalent: 

1. The function u, considered as a measurable function on N , is a finite 
energy solution to the Helmholtz equation 137]) with inhomogeneity f 
and Dirichlet data h in the sense of Definition \3.4\ 

2. We have 

{Ag + k'^)ui = fi on Ml, u\qm = h:=ho F (56) 

and 

{Ag + e)u2 = f2 on M2. (57) 

Proof As in the proof of Theorem I3.5[ we first prove that Helmholtz on M 
implies Helmholtz on N. 

Let / e L'^{M, dx) be a function such that supp (/) n (7 U dMi U SMs) = 0. 
Assume that a function u = {ui, U2) on M is a classical solution of fISBl) and 
f lFr|) . Define u = F^u and / = f oF~^ on N\T, and extend it, e.g., by setting 
it equal to zero on S. Note that then / G L'^{N, dx) is supported away of S. 
Then u E L'^{N, \'g\^^'^dx) satisfies 

{Ag + e)ui = fi = f\N, mNi, u\aN = h, (58) 

and 

{Ag + e)u2 = h = f\N2 mA^2. (59) 
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Let be the e-neighborhood of S with respect to the metric g. Let 

7i(e) be the e-neighborhood of 71 = {0} C Mi with respect to the metric 
g. Let 72(£) be the e-neighborhood of 72 = {NP} C M2 with respect to 
the metric g. Let gbnd and gbnd be the induced metrics on d'y{e) and dT,{e), 
corresp ondingly. 

Clearly, the function u satisfies conditions (120|) . (1211) . and (122|) . By Lemma 
13.71 we have that u E H^{N, \g\^^'^dx), and DgU is thus well defined. 

Using relations ©, ([6]) in Mi and ([8]) in M2, it follows from ([58]), ([59]) that 
for G Co^(Ar), 



-{Dg)udji) + fc^M^ls^l^/^ - (60) 
lim / {-f^ diudji) + {k'^u + f)^)\g\^''^dx 
lim( / + / ){-g'' V, diU^)\gbnd\^l''dS 

JdS{£)nN2 J^^E{e)nNl 

lim [ i-g'^ d,ui v,{i, o F))\gbnd\"^dS 
+ lim / i-g'^ d,U2 o F))\gbn,\^'^dS 



0. 



Indeed, both terms in the right-hand side of fl60|) tend to by the same 
arguments as the term fi'*"' diUi {tp o F))\gbnd\^^'^dS in fl28|) . Hence 

we see that (1231) is valid and thus 

{Ag + k'^)u = f inN 

in the sense of the Definition 13. 4[ 

So far, we have proven that a (classical) solution to the Helmholtz equation 
on M yields a finite energy solution to the equation on N. Next, we prove the 
converse, i.e., that the Helmholtz equation on implies Helmholtz equation 
on M. 

Assume that u satisfies Helmholtz equation flTI?]) on {N,g) in the sense of 
Definition 13. 4[ with / G L'^{N) supported away from S. In particular, u is a 
measurable function in satisfying ([T^ . flTBl) . and ffT7|) . 
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Let u = uo F and f = f o F on M \ Then we have 

{Ag + k'^)ui = fi = f\M,\^, in Ml \ 71, Ui\omi = h (61) 

and 

(A, + A;2)m2 = /2 = /|m,\7, inM2\ 72. (62) 
By conditions ([TB]) . ffTBl) . and f[T7|) . we have that 

G L\Mi\-ii,\g\^'Hx), 
gi\d,u,){dku,) G Li(M, \7., |^/r/'c/x), z = 1,2. 

Thus Mj G H^{Mi \ 7j, rfx). As before, we see that then 

(Ag + /c^)Mi = /i in Ml, Ui|aMi = /i, (63) 
{Ag + k'')u2 = f2 inM2, 

where is extended to have the value at 7^ and Ui are smooth near 7j. 
Since m satisfies ( !23|l . we see that for ip G C^{N), 



Jn 

= lim [ {-g'^diudjip + {k'^u + f)ip)\g\^^^dx 
= hm ( / +/ ) (-5'^-' <9iu|aE(£) i^i^)|5'6ndr/^c?^(a;) 

'^^^ \JdT,(e)nN2 JdT,(e)nNi/ 



'aS(e)nAf2 ^aS(e)nAfi 

= hm / {-g'^ diUiUjilj)\gbnd\^/'^dS{x) 

+ hm / {-g'J diU2iyjip)\gbnd\^^'^dS{x) 
^972(6) 

= 0, 

where ip = ip o F. Here as in the proof of Proposition 13.51 we use the fact 
that Ml is smooth function implying that diUi is bounded. 

Thus, we have shown that the function m is a classical solution on M of 

(Ag + fc2)ui = /i in Ml, u^\9M, = h (64) 
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and 

{Ag + e)u2 = f2 inMs. (65) 

This proves the claim. □ 

Next we prove a result that is not necessary for the proof but gives, in the 
case of the double coating, an alternative treatment of the distribution DgU, 
simpler than before. 

Lemma 3.9 In the double coating construction, the term 

\g\^''^g'^d.iUeV\N,dx), (66) 

appearing in Definition \ 3.4\ as DgU is well-defined as a sum of products of 
Sobolev distributions and Lipschitz functions. 

Proof. The problem we need to consider is here is that L'^{N, \'g\^^'^dx) 
contains functions that are not locally integrable with respect to measure dx 
and thus we do not immediately see that distribution derivatives dju in N 
are well defined. We deal with this by applying condition fl22|) . To do this, 
let u = M o F : M \ 7 ^ R. Using ([20]), i^, ([22]) and changing variables in 
the integration, one sees that 

/ \g\^^^g'^(,diu)dju) dx < oo. 

As g is bounded from above and below, this implies that u G H^{M\'j, dx) C 
L^{M \ 7), dx). Furthermore, changing variables again implies that 

Jn\i: 

so that u E L^{N, det (g)^^'^dx). Now in the boundary normal coordinates 
(cu, r) near S, r(x) = dist]R3(x, S), we have 

T--2|^|l/2 g [ci,C2], Ci,C2 > 0, 

and thus 

/ \u\dx = / \u\t{x)^^^t{x)'^^^ dx 

J N J N 

< \\uT^^^\\L(i{N,dx)\\T{x)'^^^\\L6/5(^N^ax) 

< ||^^||L6(Ar,T2dx')l|T(2;) '^^^ \\L^{N,dx) 

< C\\u\\L(^{N,\g\^/^dx) < 
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cf. the discussion at the end of §§3.2.3. A similar computation shows that 
u G LP{N, dx) for some p > 1, and thus djU G W~^''p{N, dx). As is shown at 
the end of §2 that 

multiphcation by \g\^^'^'g^^ maps 

H/i,p' — , H/iy and thus, by duahty, 

\g\^'^W^djU eW-^'P{N,dx), 

i.e., the distribution fl66l) is defined as a sum of products of Lipschitz functions 
and Vl^~^'^-distributions □ 



3.4 Coating with a lining: a physical surface 

In the previous sections we have considered the Helmholtz equation in a 
domain iV C M^, equipped with a metric 'g that is singular at a surface S. 
Later, for Maxwell's equations, we will need to consider S as a "physical" sur- 
face, i.e., an obstacle on which we have to impose a boundary condition. To 
motivate these constructions, we consider next, for the Helmholtz equation, 
what happens when we have such a physical surface at S. More precisely, 
we consider the Helmholtz equation in the domain \ S = A^^i U A^2 where, 
on the both sides of the boundary of S, that is, on S+ = dNi \ dN and on 
S_ = dN2-, we impose a degenerate boundary condition of Neumann type. 
In physical terms, this corresponds to having a material, sound hard surface 
located at S, separating space into two open components, Ni and A'2. Al- 
though we will not need this, it can in fact be shown that m is a solution in 
the sense of Def . 13.101 iff it is in the sense of Def . 13. 4[ 

3.4.1 Weak solutions for the double coating with Neumann boun- 
dary conditions 

In the following, we consider a double coating (M, A^, F, S, g). Suppose that 
A; > and /G L'^{N, \g\^''^dx). 
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Definition 3.10 We say that u is a finite energy solution of the boundary 
value problem with degenerate Neumann boundary conditions at S, 

{Ag + k'^)u = f znA^\S, (68) 

u\dN = h (69) 

\9\'^^d,u\^^ = 0, \g\^^''d,u\^_ = 0, (70) 

if u is a measurable function in N \ Ti such that 

ueL\N\E,\g\^/^dx); (71) 

djUeHUN\J:,dxy, (72) 

/ \g\^^^ g'^ diudjU dx < oo; (73) 
Jn\t, 

{Ag + k'^)u = f in some neighborhood of dN, (74) 
u\dN = h; 



and finally, 



for all 



[ (-g'^ diudji) + (A;2 - f)uij\ \g\^''^dx = (75) 



^i{x), X G A^i, 

V^2(x), X e N2, 



with ipi G C°°(A^i) vanishing near the exterior boundary dN = dNi \ S and 

^2 e C°°(iV2). 

Invisibility for tlie double coating with a physical surface at S, with respect 
to the Helmholtz equation at all frequencies, is a consequence of the following 
analogue of Theorem 13.81 : 

Theorem 3.11 Let u = (mi,M2) : M \ 'J ^ R and u : N \ ^ R be 

measurable functions such that u = u o F . Let f = (/i, /2) : M \ 7 ^ M 
and f : N \ Ti —>■ R be L"^ functions supported away from S and 7 such that 
f = foF, andh:dN^R,h: dM^ ^Rbe such that h = hoFi.. 

Then the following are equivalent: 
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1. The function u, considered as a measurable function on N \ is a 
finite energy solution of with Neumann boundary conditions at E 
and inhomogeneity f in the sense of Definition \3.1di 



2. The function u satisfies 



{Ag + k'^)ui = fi on Ml, u\dMi = h := ho F 



(76) 



and 



{Ag + P)u2 = f2 on M2. 



(77) 



Proof. The proof is identical to that of Theorem I3.8[ 



□ 



Remark. Let ^be a singular metric on corresponding to a double coating. 
The implication of Theorems 13.81 and 13.111 is that the solutions m in \ E 
coincide in the following cases: 

1. We have the metric g on N, singular at the virtual surface E. 

2. We have the metric g on N \J] and a sound hard physical surface at 
E, in the sense of Definition 13.101 

Similar results can be proven when the metric in corresponds to a single 
coating. 

4 Maxwell's equations 
4.1 Geometry and definitions 

Let us start with a general Riemannian manifold (M, g), possibly with a non- 
empty boundary, and consider how to define Maxwell's equations on M. We 
follow the treatment in [KLSj . using, however, slightly different notation. 

Using the metric g, we define a permittivity and permeability by setting 
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on M. 
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Although defined with respect to local coordinates, e and /i are in fact in- 
variantly defined, and transform as a product of a (+1)— density and a con- 
travariant symmetric two-tensor. 

Remark. In with the Euclidean metric gjk = Sjk, we have e^'^ — ji^^ — 
6^^. If we would like to define a generalization of isotropic media on a general 
Riemannian manifold, it would be as 

t,^'^^a{xM/'g^\ 

on M, where a{x) is a positive scalar function. However, in the following we 
assume for simplicity that a — 1. 

In the following we consider the electric and magnetic fields, E and H, as 
differential 1-forms, given in some local coordinates by 

E = Ejdx^, H = Hjdx\ 

and J, the internal current, as a 2-form. 

Now consider the time harmonic Maxwell's equations on (M, g) at frequency k. 
They can be written invariantly as 

dE = ik *g H, dH = -ik*gE + J (78) 

where *g : C°°(1]^M) — > C^iSf-m) denotes the Hodge-operator on 
j-forms, < J < 3, given on 1-forms by 

*g{Ejdx^) = ^\g\'/^g^'EjSipgdxPAdxi (79) 

= ^s^^Ej sipgdx^ A dx"^ 

where si^q is the Levi-Civita permutation symbol, si^q = 1 if {l,p,q) even 
permutation of (1, 2, 3), sipg — —1 if {l,p,q) odd permutation of (1, 2, 3), and 
zero otherwise. Thus 

*g{Ejdx^) = (e^^Ej) dx^ A dx'^ - (e^^Ej) dx^ A dx^ + {e^^Ej) dx^ A dx'^ . 

Next, we want to write these equations in arbitrary coordinates so that they 
resemble the traditional Maxwell equations. The idea is that we want to have 
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expressions which speciahze, in the case of the Euchdean metric on M^, to 
expressions involving curl and the matrices e^'' and n^''. To write equations 
in such a form, let us introduce, for H = Hjdx\ the notation 



The exterior derivative 



(curl if)' = s^^'—Hg 



OH ■ 

d{Hjdx^) = -^-^ dx^ A dx^ 

may then be written as 

dH = ^{cni\ HY sipgdx^ A dx'^. (80) 

Combining fl79|) and flHOj) we see that Maxwell equations flTHl) can be written 

as 

(curlE)' = ikfi^^Hj, 
{cnilHy = -ike^^Ej + jK 

Below, we denote also 

(V X Ey = {cnil Ey, 

and usually denote the standard volume element of by dVo{x). 

There are many boundary conditions that makes the boundary value problem 
for Maxwell's equations on a domain well posed. For example: 



• Electric boundary condition: 

u X E\qm = 0, 

where u is the Euclidean normal vector of dM. Physically this corre- 
sponds to lining the boundary with a perfectly conducting material. 

• Magnetic boundary condition: 

z/ X H\dM = 0, 

where u is the Euclidean normal vector of dM. In other words, the 
tangential components of the magnetic field vanish. 
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Soft and hard surface (SHS) boundary condition [HLSt [Kill IKi2l [Li] : 



(■E\gM = and (■H\gM = 



where C = is a tangential vector field on dM, that is, ^ x = 0. In 
other words, the part of the tangential component of the electric field 
E that is parallel to ( vanishes, and the same is true for the magnetic 
field H. This can be physically realized by having a surface with thin 
parallel gratings [HLHl IKITl lKl2l [Li] . 

4.2 Definition of solutions of Maxwell equations 

Assume that k G ]R\ {0}. We will define finite energy solutions for Maxwell's 
equations in the same way for both the single and double coatings. 

Let (M, A^, F, 7, S, g) be either a single or double coating construction, as in 
[121 denoting as usual g = F^g on N \T,. On M and \ S, we then define 
permittivity and permeability tensors by setting 



Let J be a smooth internal current 2-form on M that is supported away from 
dM. 

4.3 Finite energy solutions for single 
and double coatings 

The definition of finite energy solution is the same for both coatings. On M, 
the parameters e and fi are bounded from below and above, so Maxwell's 
equations. 



are defined in the sense of distributions in the usual way. Here, u denotes the 
Euclidian unit normal vector of dM and i? (•,-,■) is a boundary value operator 
corresponding to the boundary conditions of interest, e.g., -R(z/, E, H) = uxE 
for the electric boundary condition. 




g\'/'g^\ onAr\S. 



V X E = ikfi{x)H, V X H 
R{iy,E,H)\9M = b 



ike{x)E + J in M, 



(81) 
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If J is smooth, Maxwell's equations imply that E,H E C°°{M). 

Next, we consider Maxwell's equations on N. Let J be a smooth 2-form on 
that is supported away form dN U S. 

Definition 4.1 Let [M, N, F,'~f,Ti, g) be either a single or double coating. 
We say that {E,H) is a finite energy solution to Maxwell's equations on N, 

V X E = ikJl{x)H, V X H = -ike{x)E + J on N, (82) 

if E, H , D := e E and B := ]1H are forms in N with L^{N, dx) - coefficients 
satisfying 

ll^ll V,l?|Wo(.)) = / e''E,EkdVo{x)<oo, (83) 
ll^lli^Wi^li/WoW) = / fi''H,HkdVo{x)<oo; (84) 

J N 

{E,H) is a classical solution of Maxwell's equations on a neighborhood U C 
iV ofdN: 

V xE = ikJl{x)H, V X H = -ike{x)E + J in U, 
Riu,E,H)\aN = b; 

and finally, 



{{V xh)-E- ikh ■ Jl{x)H) dV^^x) = 0, 

N 

[ ((V X e) • ^ + e ■ {ike{x)E - J)) dVo{x) = 
for alle,h E C^{n^N). 

Here, C^IQ-^N) denotes smooth 1-forms on whose supports do not inter- 
sect dN, and the inner product "■" denotes the Euclidean inner product. 

Remark. The fact that E, H are solutions of (!82l) in the sense of Def. 14.11 
implies that they are distributional solutionsin the usual sense. Thus they 
also satisfy the divergence equations, 

V-eE = \v-J, V-J1H = 0, (85) 

ik 

in the sense of distributions. 
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5 Full wave invisibility for the double coating 



In this section, {M, N, F, g) denotes a double coating construction. In- 
visibility for active devices enclosed in the double coating, with respect to 
Maxwell's equations at all frequencies, is a consequence of: 

Theorem 5.1 Let E and H be 1-forms with measurable coefficients on M\^ 

and E and H he 1-forms with measurable coefficients on N \ 'E such that 
E = F*E, H = F*H. Let J and J be 2-forms with smooth coefficients on 
M \'-f and \ S that are supported away from 7 and S. 

Then the following are equivalent: 

1. The 1-forms E and H on N form a finite energy solution of Maxwell's 
equations 

V X E = ikJi{x)H, V X H = -ike{x)E + J on N, (86) 
R{i^,E,H)\9N = b. 

2. The 1-forms E and H on M satisfy Maxwell's equations 

V X E = ikfi{x)H, V X H = -iks{x)E + J on Mi, 
R{iy,E,H)\aj, = b 

and 

V X E = ikii{x)H, V X H = -iks{x)E + J on M2. 

Proof. First we prove that Maxwell's equations on M imply Maxwell equa- 
tions on N 

Assume now that the 1-forms E and H are classical solutions of Maxwell's 
equations on M = Mi U M2, 

V X E^ ikn{x)H, V xH ^ -ike{x)E + J on M = Mi U M2, 
R{v,E,H)\Qj,^b. (87) 

Since J vanishes near 7, ellipticity implies that E and H are smooth near 7. 
Define on AT \ E the forms E = {F-^fE, H = {F-^fH, and J = {F-^)*J. 



41 



Then E satisfies the Maxwell's equations on \ E, 

V X E = ik'jl{x)H, V X H = -ike{x)E + J on iV \ E, (88) 

Again, let S(t) be the t-neighborhood of S with respect to the metric ]) and 
7(t) the t-neighborhood of 7 with respect to g. Let If : d'~f{t) — M be the 
identity embedding. Denote by u be the unit normal vector of 9S(t) and 
d'y{t) in Euclidean metric. 

Now, writing E = Ej{x)dx^ on M, we see using the transformation rule for 
differential 1-forms that the form E = {F~^)*E is in local coordinates is 

E = Ej{x)dx^ = {DF-y.{x) Ek{F^\x))dx\ x G \ S, 

and, using Ft = F o : d-y(t) <9S(t), we have 

r{Ej{x)dx^) = {DF-^ f.{x) Ek{F-\x)) dx\ x = F{x) (89) 

Let us now do computations in the Euclidean coordinates. In the Euclidean 
metric ge, the matrix DFf^ satisfies 

\\DFf^\\(^TdT:(t),ge)^iTd-f(t),ge) < Ct, (90) 

and since E is smooth near 7 we see 

W X Eiy)\u3 < Ct, yedEit). 
Thus using §8^ we see that for h e C^{Q^N) 

[ {{V xh)-E -ikh-]l{x)H)dVo{x) (91) 

JN 

= lim / {{V xh)-E- ikh ■ Jl{x)H) dVoix) 
= -lim/ {u X E)-hdS{x) = 0. 

Thus, we have shown that 

V X E = ikJi{x)H in N (92) 
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in the sense of Definition I4.1[ Similarly, we see that 

V X ^ = -ike{x)E + J in iV (93) 
in the same finite energy sense. 

Next we show that Maxwell's equations on implies Maxwell's equations 
on M. Let U C M be a bounded neighborhood of 7 and W = F{U \ 7) U S 
be a neighborhood of S such that supp (J) fl = 0. 

Assume that E and H form a finite energy solution of Maxwell's equations 
(!86l) on (A^, g) in finite energy sense with a source J supported away from S, 
implying in particular that 

P'^EjEk e L\W, dx), p'^HjHk G L\W, dx). 
Define E = F*E, H = F*H and J = F* J on M \ 7. We have 

V X E = ikn{x)H, V X H = -ike{x)E + J in M \ 7 

and 

e^'E.E'k eL\U\ 7, dV^{x)), fi^'H.lh G L\U \ 7, dVo{x)). 

As e and /i on M are bounded from above and below, these imply that 

V X E e L\U\-f, dVo{x)), V X H e L\U\-f, dVo{x)), 
V- (eE) = 0, V- {12H) = in t/ \ 7. 

Let E'^, G L'^{U, dVo{x)) be measurable extensions of E and H to 7. Then 

V X E^ - ikfi{x)H^ = inf/\7, 

V X - ikfi{x)H' G if~^(f/, 

V X + ike{x)E^ = in t/ \ 7, 
Vx + ike{x)E' G //"^(f/, rf^ola;))- 

Since 7 is a subset with (Hausdorff ) dimension 1 of the 3-dimensional domain 
U, it has zero capacitance. Thus, the Lipschitz functions on U that vanish 
on 7 are dense in H\U), see |KKMl Thm 4.8 and remark 4.2(4)], or jXF| 
Thm. 3.28]. Thus there are no non-zero distributions in H~^{U) supported 
on 7. Hence we see that 

W X E^- ik^i{x)H^ = 0, V X + ike{x)E^ = in f/. 
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This also implies that 



V- {eE') = 0, V- ifiH^) = inU. 
These imply that E'^ and are in smooth in U. 

Summarizing, E and H have unique continuous extensions to 7, and the 
extensions are classical solutions to Maxwell's equations. 



6 Cauchy data for the single coating must 
vanish 

In this section {M, N, F,'y,'E, g) denotes a single coating construction. The 
following gives the counterpart for Maxwell's equations of the hidden Neu- 
mann boundary condition on dM2 that appeared for the Helmholtz equation. 

Theorem 6.1 Let E and H be l-forms with measurable coefficients on M\'y 
and E and H be l-forms with measurable coefficients on N \ 'E such that 
E = F*E, H = F*H. Let J and J be 2-forms with smooth coefficients on 
M \ 7 and N \ that are supported away from 7 and S. 

Then the following are equivalent: 

1. The l-forms E and H on N satisfy Maxwell's equations 

V X E = ikJl{x)H, V X H = -ike{x)E + J on N, (94) 
z/ X E\9N = f 



in the sense of Definition\4-. 1 



2. The forms E and H satisfy Maxwell's equations on M, 

V X E = ikn{x)H, V X H = -ike{x)E + J on Mi, (95) 

V X E\qMi = f 

and 

V X E = ikfi{x)H, V X H = -ike{x)E + J on M2 (96) 
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with Cauchy data 



V X = H\aM^ = (97) 

that satisfies b^ = b^ = Q. 
Moreover, if E and H solve 1^3]) . / fP^) . anc? wi/i non-zero b^ orb^, then 



the fields E and H are not solutions of Maxwell equations on N in the sense 



of Definition \4.1 



Proof. Assume first that the 1-forms E and H are classical solutions of 
Maxwell's equations in M. Moreover, assume that both E and H satisfy 
homogeneous boundary condition 

P X E\qm2 =0, Z/ X H\gM2 = 0, (98) 

that is, for the field in M2 the Cauchy data on dM2 vanishes. (Here, z/ again 
denotes the Euclidean unit normal of these surfaces.) 

Again,^define on iV\S forms E = {F~^yE, H = {F-^yH, and J= {F-^yj. 
Then E satisfies Maxwell's equations on \ S, 

V X E = ikJl{x)H, V X H = -ike{x)E + J in A^ \ E, (99) 

Again, let S(t) be the t-neighborhood of S in the ^f-metric and 7(t) be the 
t-neighborhood of 7 in the (yf-metric. 

Arguing as in fl90|) and below, we see that 

|z/ X E{y)\R3 < Ct, y E aS(t) n N2. (100) 

Recall that T,i{e) = Ni n T,{e). Then, using we see that for h e 
! {{V xh)-E -ikh-Ji{x)H)dVo{x) (101) 

J N 

= lim / {(V xh)-E -ikh-Jl{x)H)dVo{x) 
= -lim [ (u X E)-hdS{x) - [ {v xE)-hdS{x) =Q 



asi(t) JdM2 
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where we used f lTUUD and CT. 
Thus, we have shown that 

V X E = ikJl{x)H on (102) 
in the sense of Definition I4.1[ Similarly, we see that 

V X H = -ike{x)E + J on N, (103) 
also in the sense of Definition 14.11 

Next we show that Maxwell's equations on imply Maxwell's equations on 
M. 

Assume that E and H form a finite energy solution of Maxwell's equations 
dMD on {N,g). Again, define on M \ 7 forms E = E*E, H = E*H, and 
J = F*J. 

As before, we see that E and H satisfy Maxwell's equations on Mi\7i and the 
E and H are in L^(Mi, dVo{x)). Using the removable of singularity arguments 
as in the case of double coating, we see that E and H have extensions E'^ 
and H'^ in Mi that are classical solutions of 

V X E^ - iki2{x)H^ = on Ml, (104) 

V xH^ + ike{x)E^ = J on Mi. (105) 

Note that (I104p implies that, for the original field E, 

lim / {u X E)-hdS{x) =\im I {u x E) ■ hdS{x) = (106) 

JdJ:{t)nNi Jd'y{t)nMi 

where h = F*h. 

Moreover, Maxwell's equations hold in the interior of M2: 

V X E -ikfi{x)H = 0, W X H + ike{x)E = J on M2. 

Let us start to analyze what the validity of the equation V xE — ik'jl{x)H = 
on A^ in the sense of Definition 14. II implies about the boundary values on dM2. 
Using (OroSD, we see that for h G C^{n^N) 

= {{V xh)-E -ikh-Jl{x)H)dVo{x) (107) 
Jn 
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lim 



{(V xh)-E- ikh ■ Ji{x)H) dVo{x) 



(7V\Si{t) 



lim 







9Ei(t) 

vxE) -hdSix). 



uxE) ■hdS{x) + I (uxE)- hdS{x) 

dN2 



(108) 



dN2 



This shows u x E\qm2 = 0. Similarly, the equation V x H + ike{x)E = J 
holding on in the finite energy sense implies that u x H\qm2 = 0- ^ 

Assume that E and H satisfy the time-harmonic Maxwell's equations on 
M2 C such that the Cauchy data (z/ x E\qm2,^ x H\qm2) vanishes. By 
continuing E and H by zero to \ M2 we obtain solutions of Maxwell's 
equation in R^. Thus J must be a current for which there exist solutions of 
Maxwell's equations in both satisfying the Sommerfeld radiation condi- 
tion and vanishing outside A^2- Such currents are nowhere dense in L^(iV2), as 
then the fields E and H corresponding to J satisfy the Sommerfeld radiation 
condition and, using Stokes' theorem, we see that the source J is orthogo- 
nal to all (vector-valued) Green's functions Ge{- ,y,k;a) with |/ e \ M2 
and a G M^. Here, the Green's function (Ge(- , y, k; a), Gh{- , y, k; a)) satisfies 
Maxwell's equations in with current aSy and the Sommerfeld radiation 
condition. 

We thus conclude that finite energy solutions to Maxwell's equations on 
with the single coating exist only if the Cauchy data x E\qm2, ^ >^ H\qm2) 
vanishes on the inner surface of the cloaked region. Thus, finite energy 
solutions do not exist for generic sources, i.e., internal currents J, in the 
cloaked region. 



7 Cloaking an infinite cylindrical domain 

We now consider an infinite cylindrical domain, N = i?2(0,2) x M for sim- 
plicity, with the double coating. Here, i?2(0,r) C is Euclidian disc with 
center and radius r. Numerics for cloaking an infinite cylinder have been 
presented in |CPSSP] . This may also provide a picture of the cloaking that 
was physically implemented with a "sliced cylinder" geometry in [SMJCPSH] , 
although precise modelling has not been carried out. With this limitation in 
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mind, the physical interpretation of Theorems 1 7 . 1 1 and ref {single coating with 
Maxwell obstacle below is that the cloaking would be more effective with the 
insertion of a liner to implement the SHS boundary conditions which are 
necessary for the existence of finite energy solutions. 

Here, we modify the treatment from ^to the noncompact setting, blowing 
up a line and trying to obtain an infinitely long, invisible cable. 

Let 

Ml = ^2(0,2) X R, 7i = {(0,0)} X M c Ml, 
M2 = 52 X M, 72 = {NP} X M c M2 

Let M = Ml U M2, 7 = 71 U 72, 

A^i = ^2(0, 2) X R \ (^2(0, 1) X R), 
N2 = ^2(0, 1) X R, 
S = 952(0, 1) X R, 

and = ^2(0, 2) X R = iVi U iV2 U S. Let 

F = (Fi,F2) : M\7^Ar\S 

be such that 

Fi : Mi\7i-.iVi, 
F2 : M2\72^Ar2. 

are diffeomorphisms. Let X : ^2(0, 2) x R \ {(0, 0)} x R) ^ (r, 9, z) be the 
standard cylindrical coordinates on Mi. We assume that F is stretching only 
in radial direction, that is, 

X{F{X-\r,d,z))) = {F,{r),9,z). (109) 

Similarly, on M2 we have variables {r,6,z), where r = dist{x, SP) and we 
assume that F has a form analogous to f llOOp in M2. For simplicity, let gi be 
the Euclidean metric on Mi and g2 the product of standard metric on S"^ and 
standard metric of R on M2. Let g = F^g on N\'E, so that (M, A^, F, 7, S, g) 
is a double coating construction in this context. 

On M and \ S we define permittivity and permeability by setting 

^jk ^ ^jk ^ \g\V'igi\ on Ml U M2, 
~jk ^ ~jk ^ \g\y2gjk^ on \ S. 
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By finite energy solutions of Maxwell's equations on N we will mean one- 
forms E and H satisfying the conditions of Definition 4.1, where we empha- 
size the assumption that D := eE and B := JIH are in L^{N\'E, dx), making 
the integrals at the end of Definition 4.1 well defined. 

To formulate the results, we need to definethe restrictions of fields on the 
fines 7i C Ml and 72 C M2. First, assume that the 1-forms E and H on M 
are classical solutions to Maxwell's equations on M, 

V X E = iki^{x)H, in M = Ml U M2, (110) 

V X H = -ike{x)E + J, in M = Mi U M2, 
X E\qmi = /, 

where J is supported away from 7 = 71 U72. Note that then E and H are C°° 
near 7, and thus we can define the restrictions of the vertical components of 
the fields on 71 C Mi, 

C-E\^,^bl C-HU.-bt (111) 

where C := (0, 0, 1) = z := x^. 

Similarly, we can define and 62 to be the restrictions on 72 C M2, 

C-EU,^bl, C-HU.^bl (112) 
Note that 6j = bj{z) and bj = b'j{z), j = 1, 2, depend only on z. 

Theorem 7.1 Let E and H be 1-forms with measurable coefficients on M\^ 
and E and H be 1-forms with measurable coefficients on N \Ti such that 
E = F*E, H = F*H. Let J and J be 2-forms with smooth coefficients on 
M \ 7 and N \ that are supported away from 7 and respectively 

Then the following are equivalent: 

1. On N, the 1-forms E and H satisfy Maxwell's equations 

V xE^ ik'ii{x)H, V xH ^ -ike{x)E + J in N, (113) 

u X E\9N = / 
and E and H are finite energy solutions. 
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2. On M , the forms E and H are classical solutions to Maxwell's equations 
U1(J\) on M, with data 



bl = C-E\,^, bl = C-E\,„ b^,=C-H\,„ b^,=C-H\,,, (114) 
that satisfy 

bl{z) = b\{z) and b\{z) = b^{z), z eR. (115) 

Moreover, if E and H solve U10\) with restrictions (TJ^ that do not satisfy 
U15\) . then then the fields E and H are not finite energy solutions of Maxwell 
equations on N . 

Proof. First we show that the equations on M imply that the equations 
hold on A^. Assume that the forms E and H satisfy Maxwell's equations 
(11 101) in M in the classical sense, with traces (I114p that satisfy (11151) . Then 
E and H are C°° smooth near 7. 

Define 1-forms E, H and 2-form J on A^\S by E = {F'^yE, H = {F~^yH, 
and J = {{F~^yj. Then E satisfies Maxwell's equations on \ S, 

V X E = ik'jl{x)H, V X H = -ike{x)E + J mN\E. (116) 

A simple computation shows that E, H, D = eE, and B = JIH are forms on 
A^ with L^{N, dx) coefficients. Again, let E(t) be the t-neighborhood of E in 
^f-metric and 7(t) be the t-neighborhood of 7 in ^f-metric. Let It '■ d'y{t) M 
be the identity embedding. Denote by u be the unit normal vector of c?S(t) 
and i97(t) in Euclidean metric. 

Now, writing E = Ej{x)dx^ on M, we see as above using Ft = F o 1^ : 
d'y{t) — ^ 9S(t), we have in local coordinates formula (!89|) . Let us next 
do computations in the Euclidean coordinates. Using (11091) . the angular 
direction rj := de, and vertical direction ( = dz, "we see that the matrix 
DFf^{x) satisfies 

\r,-{DFi\x)ri)\^3<Ct, x G 9S(t), 

\C-{df,-\x)0\r^ = i, xeom, 

( ■ {DF-\x)7]) = 0, xed^(t), 

V {DFf \x)0 = o, xeom- 
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This implies that only angular components of E vanish on E, and we have 

\ri-E\^3<Ct, xedJ:{t), (117) 
lim C ■ ElamnN, = bp j = 1, 2, 

lim C ■ ^laswniv, = J = 1, 2 

where, for {x^,x^,x^) G S C A^, we denote 

^{x\x\x') = h]{x^), h){x\x\x^) = b']{x'), J = 1,2. 

Thus, using (ITTHD we see that for h e C^{Q^N) 

I {{V xh)-E -ikh-Jl{x)H)dVo{x) (118) 
Jn 

= lim / {{V xh)-E- ikh ■ Jl{x)H) dVo{x) 

= - lim / (z/ X E)-hdS{x), 

= - [ iux{bl-bl)C)-hdS{x) 
Jt. 

= 

where ly is the Euclidian unit normal of dN2 = E. This shows that Maxwell's 
equations are satisfied on N. Observe that if hi ^ there exists a test 
function h such that the last integral is nonzero, precluding the existence 
of a finite energy solution. Similar considerations are valid for the equation 
V X if = -ikeE + J. 

On the other hand, assume that 1-forms E and H satisfy on Maxwell's 
equations f lllSp in the finite energy sense. Then, as E and H are forms with 
L^(M)-valued coefficients that satisfy Maxwell's equations in Mi \ 71 and 
M2 \ 72, we see that they have to satisfy Maxwell's equations in Mi and M2, 
and thus they are C°°-smooth forms near 71 and 72. As and H are finite 
energy solutions on A^, the above arguments show that = 6| and fc^* = ^2- 
This finishes the proof of Theorem 17.11 □ 

Remark. If E, H, and J on M are solutions of Maxwell's equations as 
in Proposition 17.11 (1) such that conditions (114) are not satisfied, then the 
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proof of Proposition 7.1 shows that the fields E = F^,E, H = F^H, and 
J = F^:J extend to forms with L^{N, dx) coefficients that satisfy 



V X E = ikB + K, 



new 



on A^, 



(119) 



V X if 



ikD + J + J, 



new 



on 



in the sense of distributions. Here, B = JIH and D = eE are 2-forms with 
measurable coefficients and Knew = Se^E and Jnew = shSt; where (5s is a 
measure supported on E and Sg and Sh are smooth 2-forms. 

Similarly, if E, H, and J on M are solutions of Maxwell's equations as in 
Theorem 17.11 (2) with non- vanishing Cauchy data ( !97l) . we see that that E, 
H, and J on N satisfy equations (11191) with distributional sources Knew and 
Jnew defined as above. 

8 Cloaking a cylinder with the SHS boun- 
dary condition 

Next, we consider N2 as an obstacle, while the domain A^i is equipped with 
a metric corresponding to the single coating. Motivated by the conditions at 
E in the previous section, we impose the soft-and-hard boundary condition 
on the boundary of the obstacle. To this end, let us give still one more 
definition of weak solutions, appropriate for this construction. We consider 
only solutions on the set A^^i; nevertheless, we continue to denote dN = 



Definition 8.1 Let (Mi, A^i, F, 71, S, (yfi) be a single coating construction. 
We say that the 1-forms E and H are finite energy solutions of Maxwell's 
equations on Ni with the soft-and-hard (SHS) boundary conditions on E, 



ly X E\9N = /, 

if E and H are 1-forms on Ni and eE and JIH are 2-forms with measurable 
coefficients satisfying 



dNi \ S. 



V X E = ik^{x)H, V X H 
7]-E\j: = 0, r]-H\j: = 0, 



ike{x)E + J on Ni 



(120) 
(121) 



E 



L^Ni,\g\^/^dVo) 



2 




e"^ Ej Ek dVo{x) < 00, 



(122) 
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Il^ll V,|^|VWo) = / l^''H,H,dVo{x)<^- (123) 

Maxwell's equation are valid in the classical sense in a neighborhood U of 
dN: 

V xE = ikJl{x)H, W X H = -ike{x)E + J in U, 

V X E\qj^ = /; 

and finally, 

[ {(yxh)-E-ikh-]l{x)H)dVo{x) = 0, 

Jni 

[ ((V X g) ■ ^ + e ■ {ike{x)E - J)) dV^ix) = 0, 

for all e,h E C^{il^Ni) satisfying 

Ve\^ = 0, r]-h\j: = 0, (124) 
where rj = dg is the angular vector field that is tangential to E. 

We have the following invisibility result. 

In this section (Mi, A?"!, F, 71, E) is a coating configuration corresponding to 
single coating of a cylindrical obstacle -82(0, 1) x M. 

Theorem 8.2 Let E and H be 1-forms with measurable coefficients on Mi \ 
7i and E and H be 1-forms with measurable coefficients on Ni such that 
E = F*E, H = F*H. Let J and J be 2-forms with smooth coefficients on 
Ml \ 7i and Ni \ S, that are supported away from 71 and S. 

Then the following are equivalent: 

1. On Ni, the 1-forms E and H satisfy Maxwell's equations U2(J\) with 
SHS boundary conditions U21\) in the sense of Definition \8.1[ 

2. On Ml, the forms E and H are classical solutions of Maxwell's equa- 
tions, 

V X E = iki^{x)H, m Mi (125) 

V X H = -ike{x)E + J, in Mi, 

V X E\qm^ = f. 
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Proof. First, assume that the forms E and H satisfy Maxwell's equations 
(11251) in Ml. Then E satisfies identities (11171) . Considerations similar to 
those yielding formula (IllSp imply that 

[ {{V xh)-E ~ikh-Ji{x)H)dVo{x) (126) 
= lim / {{V xh)-E- ikh ■ Ji{x)H) dVa{x) 

JNi\T.{t) 

= - lim / (z/ X E) -hdSix), 

= -lim [ {u X {{ri-E)ri + {C-E)C)-hdS{x), 

JdT.{t)r\Ni 

= 

for atest function h satisfying (I124p . 

Similar analysis for H shows that 1-forms E and H satisfy Maxwell's equa- 
tions with SHS boundary conditions in the sense of Definition I8.1[ 

Next, we show that equations on A^^i imply equations on Mi. Assume that 
1-forms E and H satisfy Maxwell's equations with SHS boundary conditions, 
and internal current J, in the sense of Definition 18. 1[ Then E and H are 
classical solutions of Maxwell's equation in Mi \ 71. Let U C Mi be a 
neighborhood of 71 and W = F{U \ 71) U S be a neighborhood of S in A^i 
such that supp (J) (IW = ^. Then we have 

P'EjEk e L\W, dVoix)), Jl^'HjHk G L\W, dVoix)). 

Define E = F*E, H = F*H and J = F*J on Mi \ 71. Again, we see that E, 
H, and J satisfy Maxwell's equations on f/ \ 7, and as above we see that E 
and H have measurable extensions on 7, E'^, G L'^iU^ dVQ^x)) , such that 
V X — ikjj{x)H'^ and V x H'^ + ike{x)E'^ are distributions in H~^{U, dVo) 
supported on 71. As before, we see obtain 

V X E'^ - ik^i{x)H^ = 0, V xH^ + ike{x)E'' = in [/. 

This shows that E and H are classical solutions of Maxwell's equations on 
Ml. □ 

Similar analysis can be done in the case when we have a physical surface 
S = 5*^ X M dividing into two regions, having the SHS boundary conditions 
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on both sides, and we define the material parameters according to double 
coating construction, i.e., on both sides of the surface. 

9 Appendix: Single and double coating for 
arbitrary domains and metrics 

The constructions of §2 and the results that follow easily extend to general 
domains and metrics. Let us assume that i7 C now is an arbitrary domain 
with smooth boundary, equipped with an arbitrary smooth Riemannian me- 
tric, g = gij{x). This defines the Laplace operator with, say Dirichlet 
boundary condition, cf. Remark 3.6. Choose a point O G f2 to be blown up, 
and assume that the injectivity radius 

of {^,g) at O is larger than 3a for some a > 0. Let B{0,r) denote a metric 
ball of {M,g) with center O and radius r. Introduce Riemannian normal 
coordinates in B{0, 3a) C Q : 

X = (x\ x^, x^) ^ (r, cj), r > 0, G §^ C Tq^, 

so that X = expQ^Tu). Let /(r) : [0,3a] [a, 3a] be a smooth strictly 
increasing function coinciding with r/2 + a near r = and with r for r > 2a. 

Define, in these coordinates, 

F : B{0, 3a) \ {0} ^ fi(0, 3a) \ B{0, a), (r, u) ^ (/(r), cu). 

We extend F by the identity to Q\B{0, 3a) and obtain a diffeomorphism 

Fi:n\{0} ^ Ni = n\ B{0, a). 

Consider the metric g = Fi^g in A^^i. Observe that surfaces lying at distance 
r from dB{0, a) with respect to the metric g coincide with surfaces lying at 
distance /(r) — a from dB{0, a) with respect to the metric g. Therefore, the 
directions normal to these surfaces are the same with respect to the metrics 
g and g. In particular, the direction of these normals, in the metric 'g, is 
transversal to dB{0,a). Thus, equations remain valid if we use t — a 
instead of r — L Similarly, we again have the estimate {g]^^"^ < Ci(r — a)^. 

One may also extend the double coating construction as follows. Let {D, g^,) 
be a compact Riemannian manifold without boundary, and choose a point 
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NP G D. Using Riemannian normal coordinates centered at NP, introduce, 
similar to the above, a diffeomorphism 

F2:D\ {NP} ^N2 = D\ B{NP, b), 

where we assume that 36 is smaller than injectivity radius of D. Pulling back 
the metric (7d, we get a metric qd on D\ B{NP, h) with the same properties 
near dB{NP, b) as g has near dB{0, a). 

Observe that, as we are inside the injectivity radii, dB{0, a) and dB{NP, b) 
are both diffeomorphic to with diffeomorphisms given by expo{a(^) and 
exp^p(6cc;). Thus, dB{0,a) and dB{NP,b) are diffeomorphic to each other. 
Gluing these boundaries, we obtain a smooth manifold = A^i U A^2 U 
S with a Riemannian metric singular on S which, as one approaches S, 
satisfies conditions ([5]) . This makes it possible to carry out all of the preceding 
analysis for the double coating. 

Note that if D is diffeomorphic to 5*^ (as earlier), then is diffeomorphic 
to ~ Ml. If however D has a non-trivial topology, may have topology 
different from that of fl. However, due to the full-wave invisibility, one is un- 
able to observe this change of topology from observations made at dfl. Note 
that this is in contrast to the uniqueness result that holds for C"^ Riemannian 
manifolds [LTUj . 

Similar generalizations of the single coating construction are possible when 
dD is diffeomorphic to S^. 
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